EXPONENTIAL TIME DECAY ESTIMATES FOR THE LANDAU 

EQUATION ON TORUS 



KUNG-CHIEN WU 

Abstract. We study the time decay estimates for the linearized Landau equation on 
torus when the initial perturbation is not necessarily smooth. Our result reveals the kinetic 
and fluid aspects of the equation. We design a Picard-type iteration and Mixture lemma 
for constructing the increasingly regular kinetic like waves, they are carried by transport 
equations and have exponential time decay rate. The fluid like waves are constructed as 
part of the long- wave expansion in the spectrum of the Fourier mode for the space variable 
and the time decay rate depends on the size of the domain. The Mixture lemma plays an 
important role in this paper, this lemma is parallel to Boltzmann equation but the proof 
is more challenge. 



1. Introduction and Main Result 
1.1. The Models. The generalized Landau equation reads 



(1) 



with collision operator 



d t F + £-V x F = -Q(F,F) 
F(x,0,O = F (x,O, 



Q{F,G) = V r \ [ $(e-^)TO)V € G(0-G(OV € .F(^)]de* 

Here the positive constant e is the Knudsen number, F(t, x, £) > is the spatially periodic 
distribution function for the particles at time t > 0, with spatial coordinates x G Tf, 
the 3-dimensional torus with unit size of each side and microscopic velocity (6R 3 . The 
positive semi-definite matrix $(£) has the general form 
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where -B(|£|) = |£| 7+2 , 7 > —3, is a function depending on the nature of the interaction 
between the particles, and S(£) is the 3 by 3 matrix 

m = 

The original Landau collision operator for the Coulombian interaction corresponds to the 
case 7 = — 3. 

In order to remove the parameter e from the equation, we introduce the new scaled 
variables: 

~_ 1 ~_ 1 

e e 

after dropping the tilde, the equation (JT]) becomes 

' d t F + £ • V X F = Q(F, F) , (i,U)eT^xR+xl 3 

(2) 

k F(0,x,£) = F (x,Z), 

where T^ £ means the 3-dimensional torus with size 1/e of each side. The conservation of 
mass, momentum, as well as energy, can be formulated as 

r i,£,\Z\ 2 }F(t,x,Z)dtdx=0. 



dt 



T?. ./iR 3 



As in the Boltzmann equation, it is well know that Maxwellians are steady states to the 
Landau equation (j2J). We linearized the Landau equations fl2]) around a global Maxwellian 

1 /-l£l 2 \ 

which the standard perturbation f(x,t,£) to M as 

F = M + M 1/2 / . 
It is well-know that Q(M, M) = 0, by expanding 

Q(M + M 1/2 f, M + M 1/2 f) = 2Q(M, M 1/2 /) + Q(M 1/2 f, M 1/2 f), 
dropping the nonlinear term, we can define the linearized collision Landau operator L as 

(3) Lf = 2M~ l/2 Q(M, M l/2 f) . 

Similar to Boltzmann equation, we can decompose the linear collision operator L as diffu- 
sion part and compact part. 

Lf = Af + Kf, 

the diffusion part 

A/ = M~ 1/2 V • [(^M)V(M~ 1/2 /)] , 
where the symmetric matrix is 

(4) 0(0= / $(£-£*)M(£*K*, 
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and the compact part 
where 

Z(U*) = M(£)M&)$(Z-&. 
The linearized Landau equation for f(t, x, £) now takes the form 

f d t f + £-V x f = Lf, 

(5) 

[ f(0,x,£) = I(x,£), 

here we define f(x,t,£) = i.e. G* is the solution operator (Green function) 

of the linearized Landau equation ([5]). By assuming that initially Fo(x,£,) has the same 
total mass, momentum and total energy as the Maxwellian M, we can then rewrite the 
conservation laws as 



(6) 



/ / {l,Z,\Z\ 2 }M 1/2 I(t,x,Odtdx = 



this means the initial / satisfies the zero mean conditions. 

1.2. Preliminaries. Before the presentation of the properties of the collision operator L, 
let us define some notations in this paper. For microscopic variable £, we shall use L| to 
denote the classical Hilbert space with norm 



/ 2 



^ 1/2 

2 ■ ■ ' 



the Sobolev space of functions with all its s-th partial derivatives in L 2 . will be denoted by 
The L| inner product in M 3 will be denoted by . For space variable x, we shall 
use L 2 , to denote the classical Hilbert space with norm 



l/l W) 1/2 , 



l/el 



2 



the Sobolev space of functions with all its s-th partial derivatives in L x will be denoted by 
H x . We denote the sup norm as 



L « = sup \f(x)\ . 

xeT3 /e 



For notational of simplicity, we denote {^) s = (1 + |£|) s . In this paper, we need two 
projection operators, 

(i) For any h E L|, we define the projection to the space {M 1 / 2 } as 

Po{h) — (/i, M 1/2 ) M 1/2 (£) . 
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(ii) For any vector function u G L|, we define the orthogonal projection to the vector £ as 

*<f>»=i£f- 

First, we list some important propositions about the diffusion operator A. 

Proposition 1. For any 7 > —3, we have 

(i) 

A/ = V r [0V e /] + V 5 ■ / - (£, . 

(ii) The spectrum of consists of a simple eigenvalue Ai(£) > associated with the 
eigenvector £ , and a double eigenvalue A 2 (£) > associated with eigenvector £ . Moreover, 
there are constants c\ and c 2 > smc/j t/iat asymptotically, as |£| — )■ 00 , we /lave 

Ai(o = Cl <e> 7 , A 2 (o = C2 <e> 7+2 , 

moreover 

(e,^) = Ai(oier, 

(u,t?u) = A!(oiP(e)^i 2 + a 2 (o| [*» - p(o]^| 2 , 

(hi) For any |a| > 1, 

iwoi + i^i < Ca(e) 7+2Ha| ■ 

and 



|^Ai(0| < a<£) 7+lHa| , |^A 2 (£)| < ^<£> 7+2Ha| 
(iv) (Coercivity estimate) There exists Cq, u > swc/i t/iat 

( - A(/),/> € > cojiKO^Vlll, + IK0 7/2 p(0v € /ii^ 



It is easy to see that if we want to get coercivity estimate, we need modify the decom- 
position of L. Now, we define L = A + K, where 

(7) Af = Af - voMf) , Kf = Kf + u oPo (f). 

Under above definition, we can rewrite the coercivity estimate of Landau equation as 

< - A(/),/> € > co{||<£> 7/2+1 /||^ + ll<e> 7/2 P(0V € /||i, 



(£> 7/2+1 [/3-P(0]v,/|| 2 J 



In order to estimate the Green function of the Landau equation in next section, we need 
to recall the spectrum Spec(ek), k e Z 3 , of the operator —i-irel; ■ k + L 
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Proposition 2. [3] For any 7 > —2, i/iere exists 5 > and r = r(<5) > suc/i £/ia£ 

(i) For any \ek\ > 5, 

(8) Spec(ek) C {z E C : Re(z) < -r} . 

(ii) For any \ek\ < 5, the spectrum within the region {z G C : Re(z) > —r} consisting of 
exactly five eigenvalues {o~j(ek)} 4 =0 , 

(9) Spec(ek) n{zeC: Re(z) > -r} = {(Tj{£k)} 4 j=0 , 

and corresponding eigenvectors {ej(ek)} 4 =Q , where 

3 

a,j(ek) = a j,n(i\ £ k\) n + 0(ek) A , %,2 > , 

(10) n f 

ej{ek) = ^ ej, n (i\ek\) n + 0(ek) 4 , (e,-, e fc ) ? = <5 jfc . 

n=l 

(iii) 

4 

(11) e (-i^-fc+L) t/ = Ugf + X{|£fc|<5} £ e ^(-^( e .( £ A;), f)^(ek) 

j=o 

where 1 1 11^ 1 1 £,2 = 0(l)e~ a ^ t , a(r) > 0, X{ } i> s the indicator junction. 

1.3. Main Theorem and Method of Proof. In this paper, we decompose the solution 
of the linearized Landau equation as kinetic part, fluid part and remainder part. The 
kinetic aspect of the solution is described by the diffusive transport equation: 

d t g + £ ■ V x g = Ag . 

The operator K is a smooth operator in £ variable, we use this smoothness property to 
design a Picard-type iteration for constructing the increasingly regular kinetic like waves. 

The fluid behavior is studying by constructing the Green function with the Fourier series 
in space variable x: 

f^t \ ^ 1 iwek-x+(—iTre(,-k+L)t 

S ~ / j |rjp3 I 

k€Z 1 

here the Green function is viewed as an operator for function of £. This prompts the 
analysis of the spectrum of the operator —ineC, ■ k + L. The spectrum includes five curves 
bifurcating from the origin. The origin is the multiple zero eigenvalues of L, the operator 
at k = 0. The kernel of L are the fluid variables. The fluid like waves are constructed from 
these curves near the origin. 

With the kinetic like waves and fluid like waves constructed, the rest of the solution is 
of sufficient smoothness and it has exponential time decay rate. 
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Theorem 3. For any 7 > —2, I e L| compact support in x and satisfies the zero 
mean conditions ([6]) , then the solution 

G\I = G% K I + G% F I + G^J , 

consists of the kinetic part G* K I: nonsmooth and time decay exponentially 

(12) l|G*^llx8^ = 0(l)e-°«*; 

t/ie pari G* F /: iae time decay rate depends on the size of the domain, i.e. there exist 
5, So, Co > such that 

(i) Ife>5, 

(13) ||G^/||^ = 0, 

(ii) If5o<e<5, 

(14) 11^/11^ = 0(1)6-°^, 

(iii) IfO < e < 5 Q , 

(15) = 0(l) ir ^ e -°W^; 

and i/ie smooth remainder part G' ^J: 

(16) ||G^|| HSi| = 0(l)e-°«'. 

The Mixture Lemma plays an important role to construct the kinetic like waves, it 
states that the mixture of the two operators S and K in M* (see section |2]) transports the 
regularity in the microscopic velocity £ to the regularity of the space time (x, t). This idea 
was introduced by Liu-Yu [101 dH EE2J [13] to construct the Green function of the Boltzmann 
equation. In Liu-Yu's paper, the proof of Mixture lemma relies on the exact solution of the 
damped transport equations. In this paper, we introduce a differential operator to avoid 
constructing explicit solution, this operator commutes with free transport operator and 
can transports the microscopic velocity regularity to space regularity, this idea can also 
apply to Boltzmann equation to simplify Liu-Yu's results. Actually, the Mixture Lemma 
is similar in spirit as the well-known Averaging Lemma, see [H IU Ej . These two lemmas 
have been introduced independently and used for different purposes. 

The spectrum analysis of the Landau equation was introduced by Degond-Lemou [2], 
this result parallel to the Boltzmann case done by Ellis-Pinsky [3], however, this paper 
is the first one give coercivity estimate of Landau equation (they called Poincare type 
inequality), and they also give some ideas to decompose the collection term as diffusion 
part and compact part, see (0). However, the Poincare type inequality is not enough in 
our analysis, we thank Mouhot [13] and Guo [5] give more shape coercivity estimates for 
the linearized Landau operator. 
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The analysis based on spectrum has been carried out by many authors. The exponential 
time decay rates for the Boltzmann equation with hard potentials on torus was firstly 
provided by Ukai [16]. The time-asymptotic nonlinear stability was obtained in (9j [T7] . 
Using the Nishida approach, jH] obtained the time-asymptotic equivalent of Boltzmann 
solutions and Navier-Stokes solutions. These works yield the L 2 theory, since the Fourier 
transform is isometric in L 2 . 

In this paper, we need to restrict 7 > —2, there are two reasons about this restriction, 
first, we need spectrum gap analysis, but this only holds for 7 > —2, if 7 < —2, the 
spectrum of the collision operator is not clear. The second reason is singular part of 
compact term K s , the derivative of this operator need L| norm bound, and the coefficient 
need bounded by the cut off distance, this property holds only when 7 > —2, hence we 
need this restriction. 

The pointwise description of the one-dimensional linearized Boltzmann equation with 
hard sphere was firstly provided by Liu-Yu [10J, the fluid like waves can be constructed by 
both complex and spectrum analysis, it reveals the dissipative behavior of the type of the 
Navier-Stokes equation as usually seems by the Chapman-Enskog expansion. The kinetic 
like waves can be constructed by Picard-type iteration and Mixture lemma. In this paper, 
we apply similar ideas to Landau equation on torus, we can also construct the kinetic and 
fluid like waves, which are both time decay exponentially. Moreover, the decay rate of the 
fluid like waves depends on the size of the domain. 

1.4. Plain of the Paper. The rest of the paper is organized as follows. In section (2j we 
design Picard-type iteration and Mixture lemma for constructing the increasingly regular 
kinetic like waves. In section |3l we construct the Green function of the linearized Landau 
equation on torus, then use the long wave short wave decomposition and the spectrum 
analysis to obtain time decay rate, moreover, we improve the estimate of the fluid like 
waves. 



In this section, we will use the kinetic decomposition to construct the kinetic part, and 
apply Mixture lemma to show the tail term has enough regularity in x. 

2.1. Picard iteration and Kinetic Decomposition. We rewrite the linearized Landau 
equation §5§ as 



2. Kinetic Part 



(17) 
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The operator K has some regularizing effects with respect to the microscopic variable £. 
We decompose it as K — K s + K r , K s = K Sj d and K r = K r £>: 



K.f 



X 



D 



M -i/2 (e)M -i/2 ( ^ )v? . v^z(e,e*)/(e*)^* 



K r f = Kf-KJ, 

x(r) = 1 for r G [—1, 1] , 

I supp(x) C [-2, 2] , X G C C °°(R) , X > • 
Lemma 4. For 7 > —2, w;e /icwe t/ie following smoothing properties about the operator K 

(19) \\K s f\\Hl < D^ 2 \\f\\ q , \\K r f\\ H s < C\\f\\ q . 

Proof. It is easy to see the regular part dominated by Guassian, so the estimate of regular 
part K r is obvious. For singular part K s , we introduce the Hardy-Littlewood maximum 
function: for any / G L| 



M(f) = supr~ 3 / \f&M 

r>0 J\£~£J<r 



then d^K s f can be bounded by Hardy-Littlewood maximum function M(f) 

IW(£)I < C [ |£-£*l (7+2) - 3 |/(£*M* 
J\e-e.\<D 
00 „ 

= cj2 / ie-e*i (7+2) - 3 i/(e*)Me 

n=0 ■/2-( n + 1 )£><|£-£ ) ,|<2-™D 
00 „ 

„=0 ^|?-^|<2-«D 

< CD^ +2 M(f). 
By the Hardy-Littlewood- Wiener maximal theorem, we have 

||^(/)||l|<C||/|| l| , 

and hence 

||<W|| L| < C^ +2 ||/|| L| . 
We define the operator g = §*g an d j — ®*io as follows : 

dtg + Z- v x g - A<? = 

(20) 

s(x,0,O = 0o(3,O 



□ 
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and 



(21) 



d t J + £ • V x j - Aj = K s3 



j Or, 0,0 = j'o(x,0 . 
We have the following standard estimates about diffusive transport equations: 
Lemma 5. 



(22) \n L *Ll<e- Cot , W&W^q 



x £ 

Proof. Note that 7 > —2, we have 

(23) 2 = |Tf7| L L A{9)9dxd ^ ^ ~ cM ^ 

1/e 

this proves IIS^olLaz, 2 — e ~ cot \\9o\\ 2 T2r^- For O*, we have 

x £ L, X L I ^ 

-V^ 1/e 

we get our result if we choose D sufficiently small. □ 
Now, we design a Picard type iteration, which treat the regular part K r f as the source 
term. The —1 order approximation of the linearized Landau equation f lT7|) is the damped 
transport equation 

f d t h^ + f • V x h^ - Ati-V = K„h { - X) , 

(25) 

The difference / — h^~^ satisfies the equation 

f W - h™) + C ■ V,(/ - h™) - A(/ - h™) = K(f - M- 1 )) + K r h^ , 

(26) 

[ (/-/ i (- 1 ))(x,o,e) = o, 

this means the zero order approximation h^ ' is 

f <%/i<°) + e • - A/i(°) = AT./iC- 1 ) , 

(27) 

[ /i(°)(x,0,O = 0, 

for this process, we can define the kfi 1 order approximation /i^, k > 1 

f Sfc/iW + f • V x /i (fc) - A/i( fc ) = KhS k ~^ , 

(28) 

[ /i( fc )(x,0,O =0. 
This means the solution / can be rewritten as 

/ = + fc(o) + h 0) + . . . . 
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The following lemma gives the L X L^ estimate of h^\ 
Lemma 6. For j > — 1, we have 

(29) II^IIl^I <P +1 e-^\\I\\ Liq . 

Proof. This lemma can be proven by induction. The case j = — 1 immediately from Lemma 
[5j Now, we prove the case j = 0. By Duhamel principle, 

(30) 



It is easy to see that 



h<®(x,t,£) = / §*" Sl i ; r r .O Sl /(-,s 1 )ds 1 . 
Jo 

\\h^ LlL2 <0(l)te-^\\I\\ LlL ,. 



Assume it holds for j, then by the compactness property of K, we have 



\\h« 



E> t - s (Kh {j) )(-,s)ds 



Liqds 



Now, we can define the kinetic decomposition as 



□ 



(31) G* J = hU) +7Z = G U J + ^ , 

i=-i 

then 7?. satisfies the equation 

( d t TZ + i- V X K = LK + Kh^ , 

(32) 

( K(x, 0,0=0. 

By (|32|) . the regularity of 1Z in x is equivalent to the regularity of in x, in order to 
make sure the tail term 1Z has good regularity, we need to check has regularity H 2 . 
To proceed, we need the Mixture lemma, define the Mixture operator as follow: 



(33) M$/ 



s 2j-l 



§t-^K§ Sl - S2 K§ S2 - S3 K ■ ■ ■ S 8 ^- 1 -^ KS S2j f ds 2j ■ ■ ■ ds 1 . 



o 



'0 JO 

Under this definition, we have 



(34) h {4) {x,t,0= [ M t 2 ~ S0 K r O Sa I{-,s )ds . 

Jo 

Lemma 7. [71 [TO] (Mixture Lemma) For any fo G L 2 H^, j = 1,2, we have 

(35) ||^M*/o||^<^e-^||/o|| L ^. 
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The proof of lemma [7| will be given in next subsection. We can estimate by using 
the Mixture lemma: 

(36) \\dlh^{;s)\\ m < j\a - sojV^^II^O^Jllj^dso 



< e 2 I (s- s ) gLs ||/|| L 2 L 2 



<s"e * \\I\\ L 2 Ll 

Up to now, we have the kinetic part of Landau equation, and the tail term TZ has enough 
regularity in x, the construction of fluid part and remainder part will be given in section 

El 

2.2. Proof of the Mixture Lemma. The goal of this subsection is a short and direct 
proof of the Mixture lemma. The Mixture lemma was introduced by Liu-Yu [TO], (TTJ [121 US] 
to studying the Green function of the Boltzmann equation, and the proof relies on the exact 
solution of damped transport equation. However, it may difficult to construct the exact 
solution of S*/o, in order to avoid this difficulty, we need to introduce a differential operator: 

V t = tV x + V 5 . 

It is important that T> t commutes with free transport operator: 

[D t ,d t + £-V x } = 0. 

We have the following estimates about differential operator T> t and solution operator §*. 
Now, we have the following estimates about T> t and §*. 

Lemma 8. For any f G L 2 X H^, j = 1,2, there exists rj small enough such that 

>?«3* fJI .„ .„ <r e -(<x>-no)t\ 



(37) 



ll^SVolUi^^e-^-^H/oll™, 
||2> t S*/o - §V 5 /o|| L|L | < e-^-^\\f4 m , 

k ||2^s*/o - A^Ve/olU^, < e-^-^H/oll^ 



Remark: Although the integral operator K has smoothing property, it only allows us 
differential with respect to £ once, if we calculate second order Mixture operator ^M^/o, 
the second derivative term d^K appears, so we need cancelation properties (l3"T]) q and 
( |37|K to overcome this difficulty. 

Proof. Let = 2^S*/o, 3 = 0, 1, 2, then the energy estimate for /(°) gives 



li/ (0) n^ + co jf \^y /2+1 f {0) \\l iq + iKe> 7/a p(Ov c / (0) ii^ 
+ ||<e> 7/2+1 [/3-p(o]v^||^ L? ^<ii/oiiL L? . 
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We study the operator T> t componentwise, on the other hand, define T> it = td x i + dp 
i = 1, 2, 3. We can check the commutator 

[v it ,A]h = v ? • [(<%tf)v^] + v ? • [deW)]h- [%(£,tf£)]/> 

Denote 

5 

[V t ,A]h = Y,B l (h), B l =(B[,B l 2 ,B l 3 ), 
i=i 

then solves the equation 

ft/ (1) +e-V./W=A/W + [X>*,A]/W, 
by energy estimate of jf™, we have 



Now, we estimate B 1 , ■ ■ -,B 5 separately. For B 2 , note that |V$ • | < (£) 7 , hence 



J-r/ / bv o) ) ■ /w^de <c'ii<e) 7 V ) iii^+5 1 |i<e> 7/2 / (1) iii |i| 



For I? 3 , we have | < (£) 7+2 , this means 

— , , b 3 (/(°>)./ (1) m <^lK0 7/2+ V (0) iil 2L? + 5 2 iK0 7/2+1 / (1) iii 2L? , 




l^l/ e l A3 



and for 5 4 + 5 5 
1 




|f p / (5 4 + 5 5 )(/^)-/ (1) «( <^ll/ (0) Hi^+^||/ (1) |li. L : 



'r 3 Jt 3 , 

1/e 



Finally, for B 1 , we need to decompose the microscopic velocity £ near origin and away 
from origin. Notice that for a smooth cut-off function x{0 = 1 f° r |£| — lj x(0 = for 
|£| > 2, going back to componentwise expression Bj, we have 

/ B}(fW)fiM= [ xBl(f i0) )f { M+ [ {i-x)B}(f®)f<M- 

Jm 3 ^I5I<2 ^I£I>i 

For |£| small, we have 

x^(/W)/«de= / -v c (x/ (1) )(%^)v 5 /(°^e 

•^l£l<2 

:2 7 V ^l£l<2 7 



< c 
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For |£| large, we have 

(1 - x)Bt(fW)f {1) dt = [ -V C ((1 - x)f {1) )(d^)VJ {0) dt 
ICI>i J\z\>i 

< c\\ <e> 7/2 v € /<°> || « || <e> 7/2+1 / (1) lb + / % (i - x) Ve/Wtfv^de 

5 5 ^|€|>i 

+ / (1 - x)V c %/ (1) tfV^ ^ + / (1 - x)V £ / (1) 7?V^/C ># . 

We only need to estimate the last two integral, 



/ (i - x)VtdefW#Vif<®dt + [ (l - x)v^ (1 W € %/ (0) ^ 

J\S\>1 J\t\>i 

Jl - X ) [ai(p(Ov € %/ (1) ,p(Ov € / (0) ) + A 2 ([/ 3 - p(0]v 5 %/ (1) , [h - P(O]V«/ (0) )]de 
^(1 - X ) [Ax (p(0 v^ (1) , p(0 v^ ei / (0) ) + A 2 ([/ 3 - p(0] v c /« [/ 3 - p(0] v € <W (0) 

-/|£|>1 ■'|£|>1 

I >x (i - x)%a 2 (0 ([/ 3 - P(0]v e /«, [/ 3 - p(0] v € / (0) )de 

/ (i - x) Mo - a 2 (o] ([/ s - p(0]v c / (i) , ei ) (v ^,° a '^ de 

/ (i - x) [Ai(0 - A 2 (0] ([/s - P(O]V,/ (0) , e„) ^r^ ff • 



r lCI>i 
Hence we have 



(l-x)Bl(fM)fMdt 

ei>i 



< c(ii<e> 7/2 v,/(°)|ii i + iKe> 7/2 p(ov t /W|i^ + ii<e> 7/2+1 [/3 -p(0]v,/^ni, 
+ ^(iKe) 7/2+1 / (1) iii ? + \m J/2 no^f {1 % + n<e> 7/2+1 [/3 -nwtf (1) \\l{) ■ 

Summing up above estimates for B\, B\ and B\, we have 

< ^(iKe> 7/2 p(e)v,/(°)|i^ + ii<e> 7/2+1 [/a -p(0]v^°>u^) 
+ 3 r (n(e> 7/2+1 / (1) iii i + ii<e> 7/2 p(ov € /«ii^ L| + ii<o 7/2+1 [/a -m]^f {1) \\l iq 
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where 6 depend on 64,65. Choose 61,62, 63 and 6 small enough, then there exists 770 small 
enough such taht 

(38) < "(co - Vo) f \\f (l) \\%qds + H/olli^ + IIVt/oH^ , 

X £ X £ X £ X £ 

and hence we obtain 

For second order estimate, we have 

(39) cV {2) + £ • V x f^ = A/ (2) + 2 [A, A]/W + [P t , [D t , A]]/(°) , 
the second order commutator \D t , [D t , A]]/i = {^4ij(^)}f J= i is a matrix, where 

- u (M 1/2 , ^i^h)^M 1/2 + v^M^^^M 1 ' 2 ) , 

similar argument can get our result and hence we omit the detail, this proves the inequality 
(|37|li. For dSZDa and (j37j) s , we can define g = V 5 / and g& = V^go, then gW and #W 
satisfy the following equation respectively 

% (0) + £ • V,<? (0) = Ay® , 



(40) 
and 
(41) 



/»(x,0,O = V s / , 



^)(x,0,O = V|/ . 

We can define u w = - j = 1,2. Then u« = P t S*/ - S'V^/o and w (2) 

I> t 2 §*/o - X>iS*V s / - For (|37|)s. note that solves the equation 

cV 1} + £ • V^ 1 ) = AuM + [V t , A]/(°) , 



(42) 

the energy estimate gives 



u 



(1) (^o,o = o, 



\LIL\ » 



ll« (1) ll£»za <-(<*>-»») /V (1) lli^^+ll/ol 

this complete the proof of fl37|) -). For (j3"7|U . u^ 2 ^ solves the equation 

r ftti (2) + e . v^ 2 ) - A M ( 2 ) = 2[V t , A]/W + [X> t , [A, A]] /(«» - [A, A]^°) , 

(43) 

[ u( 2 )(x,0,O = 0, 

similar argument can get our result. □ 
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Proof of the Mixture Lemma. For j = 1, we can write down d^M^fo as follows: 

d 1 x M[f (x,O= [ I 1 d x ^- Sl KS s ^K^f ds 2 d Sl 
Jo Jo 

Jo Jo Jm 6 

Note that [<9 X ,§*] = 0, [9 X ,W] = 0, we can change the order of (d x ,E> t ) and (d x ,W). In 
order to get time integrability, we can rewrite dlWl\f Q (x, £) as 

diM\f (x,o= f [ L I £i ^^ si ^(e,ei)^s si - s w(ei,6)s s Vo«i^ 2 ^ 1 

Jo Jo Jm 6 s i 

Jo Jo Jm. 6 s i 
Using the fact £V X = T> t — V^, we have 

d i x M t 1 f (x 1 o= [ r [ -§*- s w(e,ei)(^ sl - S2 - v 6 )§ si - s w(ei,6)s S2 /o«i«si 

Jo Jo Jm 6 s i 

+ [ f 1 f -S'- s W(e,ei)§ Sl " s W(6,6)(^ 2 -V 6 )§ S2 /orf«i^ 2 rf Sl . 
Jo Jo Jm 6 s i 

By (|37|) i and integration by parts, we have 

Wd'MfoWLlLl < e-^(||/o||L^ + \\dlfo\\ m )( f f 1 -ds 2 d Sl 
1 Vo JO s i 

= te- 2 -¥{\\fo\\ Ll q + \\dlf \\ Liq ). 
This proves the case j = 1. For j = 2, we can write down c^M^/o as: 
^M f 2 / (x,0 




If 7R6X2 



&- si w(t, 6)s si - s w(ei, 6)s S2 " s w(6, e 3 )s S3 - s w(e 3 , £ 4 )§ S4 /o 



dEdS 



where 

dS = dsi ■ ■ ■ ds^, dE = d£x ■ ■ ■ d£ 4 , T = [0, t] x [0, Si] x • • • x [0, s 3 ] . 
In order to get time integrability, we need decompose sis 3 as 

S1S3 = [Ol - S 2 ) + (S 2 - S3) + (S3 - s 4 ) + s 4 ][(s 3 - s 4 ) + s 4 ] . 

We then have 

d 2 Mf (x,O= [ [ —(Jifo + J2fo)d~dS, 

where Ji collects all the terms that each W differential with respect to £ at most once, and 
J 2 collects all the terms that one of W differential with respect to £ twice. More precisely, 

j 2fo = §*-w(e,6)s si - s w(ei,6)[^ 2 - S3 s S2 ^ 3 ] [v i3 w^ 2 ^ 3 )p s ^w(^uP S4 fo 
+ §^w(e,6)s si - s w(6,6)s S2 - S3 [vlw(&^ 3 )p s ^w(^u)^fo 
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(44) + §*- s W(£, £i)§ sl - s W(£i, £ 2 )S S2 - S W(£ 2 , £ 3 ) [P S3 _ S4 § s 3- 4 ] ^ w ^ ^ g«4 /o 
+ § f - s W(£, a)S sl " s W(£ 1; 6)S S2 - S W(£ 2 , £ 3 )S S3 " S4 [V 2 4 W(£ 3 , &)] § S4 /o 
+ §*- s W(£, ei)§ Sl " s W(ei, 6)S S2 - S W(£ 2) £ 3 ) [P S 2 3 „ S4 S S3 " S4 ] e 4 )S S4 /o 
+ §*- s W(£, 6)S Sl - s W(ei, 6)S S2 - S W(£ 2 , £ 3 ) [P S3 -S4§ S3 - S4 ] [V^(£ 3 , £4)] S S4 /o • 

The estimate of Ji is similar to j — 1, 



For J 2 , the first four terms of (1441) can be estimated by ( 1371) 9 and the last two terms of (J4 
can be estimated by (I3TIU . then 



J 2 fodZ 



< e 3 II fnll r 2 172 
- IUU||L^ ? 



This means 



I^M*/o|| l;l? < e~r- H/olUgfli 



— dff) <t 2 e-^||/ || L 2^2, 



this completes the proof of the lemma. 



□ 



3. Fluid Part and Remainder Part 
3.1. Green function. Consider the linearized Landau equation 

d t f + Z-V x f = Lf, 
/(x,0,£) = /(x,£), 



(45) 



where / satisfies the zero mean conditions Hereafter, we will use just one index to 
denote the 3-dimensional sums with respect to the vector k = (k±, k 2 , k 3 ) G Z 3 , hence we 
set 

E= E ■ 

(fcl,fc2,fc3)GZ 3 



Consider the Fourier series of initial condition I in x 



(46) 



/vc-: 



(/)*(£) 



IT 3 



-iirek-x 



dx 
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rewrite the solution f(x,t,£) of (14"5|) as Fourier series 



17 



(47) 



fcez 



1 



■l/e\ JTf /s 

the Fourier modes of fj46|) - fH7|) satisfy the following equations 

{djk + ineZ, ■ kf k - Lf k = , 

Hence 

(49) A(*,0 = e ( -* ra «-* +L)t (/) fc (0, 

the solution of (1451) is given by 

/0M,£) = ^ e ^-+(— 5-^)t ( j )fc(0 ^ 



k& 



IT3 



iirek-(x— y)+{— iire^-k+L)t 



the Green function G* can be expressed as 



T — 

IT?, 



D iirek-x+(— ine^-k+L)t 



(50) 

fcGZ 1 J1 'l/e | 

Note that if e — >■ 0, the Green function Gq(x, £, i) becomes an integral formula 



(51) 



3 ir?-x+(-i7rg-r;+L)t 



c??7 . 



3.2. Long Wave Short Wave Decomposition. We can decompose the Green function 
fl50|) into the long wave part G* L and the short wave part G' s given respectively by 



^iirek-x+(— iire^-k+L)t 



(52) 



|efc|<<5 1 V s 1 
| £ fc|><5 ' x / £ ' 

The following long wave short wave analysis relies on spectral analysis (Proposition [2]). 
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Lemma 9. (Short wave G* s ) For any s > 0, I G H^L^, we have 

m, s i\\ Li q<e-^nnLiq, 

(53) 

Proof. Note that 

G* jS 7 = e { - i7Te ^ k+L)t e inek - x I k ^) . 

\ek\>S 

For L? X L 2 ^ estimate, we have 

\\^ s n\i = E ie ( -^- fc+L)i 4(or, 

|efc|>5 

and hence by spectrum property (Proposition [2]), we obtain 

\KsI\\' m = E l|e^ + ^/ fc (0llL<e-°^ E ||4(0II^ 



e -o(i)t 



/ E iA(oi a de<e-° (1)t ii/|i^ 



|efc|><5 

For the high order estimate, we have 

\ek\>5 

and hence 

\K S I\\ 2 m < e-° (1)t E C 1 + M| a )'||/ fc ||^ < e-W||/||^ Li . 

efc|>5 

□ 

In order to study the long wave part G* L , we need to decompose the long wave part as 
the fluid part and non-fluid part, i.e. G* L = G* F + G* L . ± , where 

i 



G ^ / = E ^e^-^e^ie^J^^iek) 
(54) I£fc|<5j =° 

&U±i = E e^n*/* • 



\ek\<S 

Lemma 10. (Long wave G* L ) For any s > 0, I G L^Lj and satisfying the zero mean 
condition we have 

\\Gl,L;J\\H°q < e-W\\I\\LlLl , 

(55) 



h, F I\\ H s L > < e-°^ H \\I\\Liq 
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Proof. For the non-fluid part, using the spectrum property (Proposition [2]), we have 

II«WII^< E (i + k^iT|iV fe (0l 2 < (i + MT E in^(Ol s 



|2 



\ek\<5 \ek\<5 

and hence 

im^Wmq z z- oil)t \\n 2 m ■ 

For the fluid part, by spectrum property (Proposition [2]) and the zero mean conditions ([6]) 
we have 

4 

\\G e , F i\\ 2 m < {i + \7rek\ 2 y Ei^^n^W'^i 3 

|efc|<<5 3=0 

^ c E E e " aj,2|fc£|2[1+0(52)1< H^lli 2 

|efc|«5 3=0 

<Ce-°^\I\\l iLr 

□ 

Remark (i) If e > 5, we do not have long wave part, i.e. G* L = 0. 

(ii) The high order estimate of short wave part requires regularity in x. This is because 
\nek\ may not be bounded. One needs the regularity of I to ensure the decay of G* s in 
time. 

(iii) In order to remove the regularity assumption in x, we need the Picard-type iteration 
for constructing the increasingly regular kinetic-like waves in section [2j 

Theorem 11. For any I G H^L^ and satisfies the zero mean conditions ((6]), we have the 
following exponential time decay estimate about the linearized Landau equation (145 p 

||G< I\\ L „ q < e~ Xst \\ I\\ m + e- A ^|| I\\ Liq . 

(i) If e > 5, then \ s = 0(1) and \ L = oo, 

(ii) Ife<5, then X s = 0(1) and \ L = 0{l)e 2 . 

3.3. Improvement of Fluid Part. In this subsection, we improve the estimate of the 
fluid part. Recall the fluid part of the Landau equation ( 154| 

4 

G IfI= E ^e^ k ^ mk - x {e j (ek),I k ) ( e j (ek). 

\ek\<8 j=0 

We have 

4 

\KfI\\Ii = J2\ E e^ ( - £fc) V^< ej (^),4>^ 

j=0 \ek\<5 
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1 



Note that 

IT 



1/*W* ' 'V lek \ <St i/2 11 JB Kl/a (0) 



as £ — > 0, this means for any a > there exists 5q > such that if e < 5q, then 



1 V" e -( £ |fc|) 2 <ao + [ e- y2 dy = C . 



IT 3 

1 1 / £ ' lefcK^t 1 / 2 



Hence 



< (nrh E e " (£|fc|)2t ) / H 7 N^ 



<o(i) 



_1 / el |efc|<5 
Co _ £ 2 t 



(1+t) 3 / 2 



Theorem 12. Assume that e < 8, I 6 L| u>rf/i compact support in x and satisfies the zero 
mean conditions ([6]), £/ie« t/iere eins£ Co,5o > suc/i that ifO<e<5o, then 

\K F n Ll <o(i)-^^e-^. 

If e — y 0, i.e. the whole space case, the pointwise estimate of the fluid part becomes 



J* J|U < O(l); 



1 

- — (1 ^tpT^ • 
This recover the whole space result in [TO] . 

3.4. Remainder Part. The remainder part can be defined as follows: 

(56) Gt >R I = TZ- Gt >F I , 

combining long wave-short wave decomposition and kinetic decomposition, we have 

4 

(57) G% R I = Gl s I + G% L .J - £ h [j) = n- Gt jF I , 

3'=-l 

hence by (j53]l, (153]) and (125]) , 

(58) ll<VLz> < t 5 e-W||J||^ 
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This means R I is in non- fluid part, hence 

(59) ^II^G^JH^ < -c\\dlQ% R I\\ Liq + \\K(dlh^)(;t)\\LiLl ■ 
We only need to estimate 

(60) f\\K{%hW){-,t- 8 )\\^ds = 0{l). 

Jo 

This immediately from (fl9l) and (J36l) . 

J* ||K(^))(-,t - < jf lU^II^II^^ 

= o(i)\\i\\ LlLl . 

Hence we have 

(61) \\G% R l\\ H2L2 <0{l)e-°V\ 

moreover, by Sobolev inequality, we have 

(62) \\Gt tR l\\ LooL2 <0(l)e-°^ 
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